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ABSTRACT 

We reanalyze precision LEP data and coupling constant unification in the minimal 
supersymmetric SU(5) model including the evolution of the gaugino masses. We derive 
general bounds on the primordial gaugino supersymmetry-breaking mass-scale m^2 in 
terms of the various input parameters. The model cannot accommodate m^2 < 1 TeV for 
values of 0:3(771 z) < 0.115, even for extreme 1 — a values of the other inputs. We emphasize 
the sensitivity of this type of calculations to the various input parameters. 
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1. Introduction 

Gauge coupling unification has always led to attractive predictions of Grand Uni- 
fied Theories, GUTs |l|. Although qualitatively successful, minimal non-supersymmetric 
GUTs predict sin 2 Qwi^z) about 8 to 5 percent lower than the present LEP value 0. 
However, minimal supersymmetric GUTs predict sin dw{mz) tantalizingly close to the 
LEP value || . Some authors have suggested that the amazingly close agreement of LEP 
data with minimal supersymmetric GUTs is evidence for supersymmetry at a scale near 
1 TeV , in accord with theoretical prejudice based on naturalness arguments. Previously, 
we developed analytically one-loop expressions including new particle thresholds, and cor- 
rected these expressions to two- loop accuracy ||. These expressions were then applied 
to to a specific model, minimal supersymmetric SU(5) ||, in an attempt to bound the 
supersymmetry-breaking scale in the form of a universal primordial gaugino mass, mi/2 

0- 

We found that the data then available favored rax/2 < 65 GeV or m^ 2 > 21 TeV at 
the one-standard-deviation (1 — a) level. However, in those papers we did not include 
the dependences of the physical gaugino mass thresholds on the gauge couplings, which 
has recently been shown to have a significant impact on the favored range of rax/2'- the 
evolution of gaugino masses (EGM) effect || . 

In this letter, we correct and continue our previous work J7|] on the derivation of 
rigorous bounds on rax/2-, incorporating this EGM effect and the latest available LEP 
data. Since the largest uncertainty in m^ 2 is now that due to a.2,{ja z), we present our 
results as bounds on m 1( / 2 as a function of a^irnz) for specific values and plausible ranges 
of the other parameters: sin 2 Qw{ m z)-, a ern (mz), m tl the Higgs superpotential mixing 
parameter \i 1 and the Higgs mass nih- Throughout this paper, sin 2 9w{f^z) and all gauge 



couplings are defined using the MS scheme. Using extreme 1 — a values of the other inputs, 
we find that mi/2 < 1 TeV can be obtained only for as(mz) > 0.115. We show how to 
modify our results for different values of the inputs, noting that improved precision will 
most likely increase the minimum value of a^(mz) compatible with mi/2 < 1 TeV. 



2. Basic Formulae and LEP data 
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The prediction for sin 9w(jnz) in minimal supersymmetric SU{h) may be written as 



sin 2 e w (m z ) = 0.2 + a e™K™z) + Qm2 g + SJlight) + 5 s (heavy) + 5 s (conv) (2.1) 

15a 3 {m z ) 



where 0.0029 corrects the analytic one-loop calculation to two-loop accuracy. The scheme 
conversion term 8 s (conv) is negligible. Because of constraints on the Higgs triplets me- 
diating proton decay, the heavy particle threshold correction 5 s (heavy) > in minimal 
supersymmetric 577(5) [0. The contribution from light particle thresholds is: 

,.(«,«) =£=£*) (-3^) + ^InC-^ll) - ^InC-^) 
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and y = (mo/mi/ 2 ) 2 , where mo is a universal primordial supersymmetry-breaking spin- 
zero mass, w was defined in 0] and the logarithms should be set to zero if the threshold is 
below mz- 

The value of sin Qw{ m z) is known from an analysis of precision data from LEP and 
elsewhere: 

sin 2 6 w (m z ) =0.2327 ±0.0007 |] (2.4) 

The value of ot, ern {mz) is also well-known: 
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"»W = l2MfT02 EM (2-5) 

However, a glance at Table 1 shows that there is a huge range in the experimental determi- 
nations of oi^[mz) ||12|| . Naively taking an equally- weighted mean of these eight different 
values and the uncertainty of the least uncertain experiment gives a^{mz) = -116 ± .005. 
However, it is disturbing that the different experimental determinations of a^mz) do not 
overlap. Except for the DELPHI value, all the LEP measurements are higher than, and 
do not even overlap the other three lower-energy measurements. It is known that higher- 
order corrections to these LEP values are important, and preliminary attempts to sum 
the dominant contributions to the perturbative higher-order effects indicate that the LEP 



results could be shifted down to values nearer the other measurements [|T3| . 

We take an agnostic approach in this paper and derive bounds on m\/ 2 for arbitrary 
values of a^{mz), and await further progress in its determination. 



3. Derivation of Bounds on m\i<i 

In this section, we use the information of the previous section to derive bounds on 
rai/2 for different values of as(mz), assuming 5 s (heavy) = 00 and ignoring the negligible 
S s (conv). 

Consider first a scenario where the winos are lighter than mz but the gluinos are 
heavier than mz- The formula (|2.1|) then gives: 



(n(^) = \x - 1- - c- 9 (3.1) 
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where 
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Now consider the case where both the gluinos and the winos are heavier than mz, in which 

case ( |2.1| ) gives: 

In^^dl) = -X + — + 7c5-8c™ (3.3) 

m z «3 



To obtain the most generous bounds from ( |3. 1|) and (|3.3| ), we would need to maximize X, 
minimize c§, and maximize c^. 

For physically-relevant values of w (those which give positive squared masses for the 
stop squarks), f(mi/2,y,w) is minimized at w = —S(cg + y)/27. With this value of w, 
f(y, w) has one extremum, a maximum, at y — (cj cj +2cj Cq — cj Cq)/(3cj — 2cf —Cq), and 
approaches -0.025 as y becomes very large. Since the values of the c's that we encounter 
satisfy Cq > 1 > cj , cj , the minimum of f(y, w) is indeed -0.025 for values of y > 0. We 
have searched numerically the region where some of the scalars contibuting to f(mi/2, y, w) 
are lighter than mz to verify that it does not take smaller values in this region. 

Taking the extreme 1 — a values sin 2 Qwij^z) = 0.2334, m t = 95 GeV, a ern (mz) = 
1/128.1 and using a naturalness bound of 1 TeV on \x and m^, gives a maximum numerical 
value for X of 193.2. Note that the extreme values of sin 2 Ow^z) and mt are correlated, 
and that the contribution of the squarks and sleptons which is represented by f(y,w) is 
negligible. 



If 5 s (heavy) > the bounds will be even tighter. 



Estimating the values of c g and c$ is more involved, since they depend on the unkown 
value of mi/2 and an iterative procedure must be used ||. At the one-loop 

__ a3(cgm 1 / 2 ) _ cx2{cwm 1/2 ) 

Ch — C w — l^-^J 
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where c g mi/2 and c^mj^ are the gluino and wino masses renormalized at m g ,m^. 

From the one-loop expression for renormalizing a coupling from mz up to its gaugino 
threshold: 

/ v a(m z ) ,, _x 

a\m g augino) ~ _ fe f ( a(m gaugino )m 1/2 ) ^ 6 '°) 

We see that, for b < 0, a(m gaugino ) increases with b. In order to minimize c g , we want 
to use the minimum value of 63 = —7 possible in the MSSM below the gluino threshold. 
Similarly, to maximize c^ we use the maximum value of 62 — — 1/3 possible below the 
wino threshold. 

Fitting the results of an analytic one-loop calculation to a numeric two-loop calcula- 
tions for central values gives |7j: 
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where the stop squarks have been taken degenerate with the other squarks. Thus, ol\j 
decreases with the thresholds. Taking upper bounds of 165 GeV on the top mass [9] and 
3 TeV on the other thresholds gives the range 

3 1 13 1 

+ - — -, r-0.7< — < — + - — - - + 1.4 (3.7) 
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for the coupling at the unification scale. Numerically, we find a slight variation of the 



solutions of (P-lj) and (3.3) over this range of au, with both values increasing with oljj. 



Therefore, we use the maximum value in fl3.1|) and the minimum value in (3.3). 

Finally, we use the central value of a ern (mz) to calculate c^ and au, neglecting the 
small variations from the uncertainty in a em (mz)- 



4. Results and Conclusions 

Fig. 1 shows on the left the triangular region excluded by the calculation of the 
previous section, which is based on the formulae of Ref. [7] with the EGM effect incor- 
porated as suggested in Ref. [8]. The upper solid line is obtained by iterating (|3.3|). 
describing a scenario where m^ > mz, and the lower solid line is obtained by iterating 
( P-lp , describing a scenario where m-Q < mz- The area between these two solid lines 
is excluded. An approximate lower bound of TO1/2 > 45 GeV from CDF gluino searches 
and an upper naturalness bound m 1 / 2 < 1 TeV are represented as dashed lines. To have 
TO1/2 < 1 TeV, a.^{mz) > 0.115 is needed. Note that the J/\I/, T and deep inelastic values 
of as(mz) are below this lower bound! If the lower bound on m 1 / 2 can be pushed to 
60 GeV, mi/2 > 1 TeV will then require a^(mz) > 0.116. 

Table 2 shows the value of X for different values of the inputs within their 1 — a range. 
Since our results use extreme 1 — a values to maximize X , X can only decrease if these 
errors are reduced , giving an even tighter lower bound on a^(mz)- For reference, we also 
show in fig. 1 the bounds from X = 189.4 corresponding to central input values as a dashed 
line, and the minimum value of X = 185.4 as a dotted line. Note these bounds are very 
sensitive to the values of the input parameters. Since the various approximations which 
have been made all loosen the bound on 0:3 (mz), the bounds would become stricter if our 
hypotheses could be tightened. Whilst the minimal supersymmetric 577(5) GUT model 
cannot be ruled out on the basis of coupling constant unification, things are becoming 
very tight. Explaining the discrepancy between the different determinations of ot^(mz) is 
crucial. If one believes the high values, or even the average values of a^(mz), the model is 
in good shape. However, if higher-order corrections [13] do actually bring the high values 
of a^(mz) down to match the low values, then coupling constant unification cannot be 
achieved in the minimal model at the 1 — a level. 

We conclude this paper with a brief explanation of our attitude towards these detailed 
estimates of the supersymmetry-breaking scale in the minimal supersymmetric SU(5) 
model. We regard these calculations as illustrative, relevant for developing the technology 
needed to calculate this scale in any explicit model, and indicating how far one can get on 
the basis of our present experimental knowledge using the theoretical tools available, but 
we do not take the minimal supersymmetric £77(5) model as seriously as do some of our 
colleagues. Even in a purely field-theoretical context, this model has to be modified if the 
Higgs doublets and triplets are to be split in a natural way, and possible difficulties with 



too-rapid proton decay are to be avoided |T^]. More fundamentally, one cannot nowadays 



be blind to the promise and successes of string theory in unifying quantum gravity with 
the other particle interactions. It is known |0| that the adjoint Higgs multiplet needed 
in the minimal supersymmetric SU(5) model cannot be obtained from string theories in 
which gauge charges appear at the Kac-Moody level 1, which is the formulation used in 
all models constructed so far. Among existing string models, the one closest in philoso- 
phy to traditional SU(5) is the flipped SU(5) model flSH , which also incorporates natural 
Higgs doublet- triplet splitting. This model offers the possibility that sin 2 Qw( m z) could 
be smaller than in minimal supersymmetric SU(5), which is not inconsistent with the 
data, particularly if the lower quoted values of as (mz) turn out to be correct. In this and 
many other string-derived models the value of the string unification scale at which all the 



gauge couplings appear to become equal is calculated [17] to be significantly higher than 
the grand unification scale calculated in minimal supersymmetric SU(5). Within the con- 
text of flipped SU(5), this reinforces the suggestion that sin Qwij^z) might be decreased. 
More generally, this scale discrepancy suggests the existence of additional particles beyond 
those in the Standard Model or its minimal GUT extension. One example, the String 



Inspired Standard Model (SISM), was proposed in ref. |fL8| . It would be desirable to carry 
out in this model and others derived from string calculations as detailed as those described 
in this and our previous papers, but such a study would take us beyond the scope of this 
note. The reader should keep in mind, though, the observation that analyses of precision 
LEP data can be used to test string ideas, as well as the GUT ideas discussed in the bulk 
of this paper. Yet more motivation, as if our experimental colleagues needed any, to carry 
on compressing the LEP error bars, in particular for 0:3 (m.z). 
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Figure Caption 

The allowed ranges of TO1/2 for different values of as(mz) are shown, in the context 
of minimal supersymmetric SU(5). The value of mi/2 for the central values of the other 
inputs is displayed as a dashed line. The minimal value of as(mz) assuming 1— a variations 
of the other inputs is shown as a solid line, and the maximal 1 — a value of as(mz) as a 
dotted line. For any given value of a^imz), values of mi/ 2 to the left of the solid line, or to 
the right of the dotted line, are accordingly disfavoured. Also shown as horizontal dashed 
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lines are a lower bound toim > 45 GeV from CDF gluino searches, and an upper bound 
777-1/2 < 1 TeV motivated by naturalness. It can be seen that mi/ 2 < 1 TeV is favoured 
only if as(mz) > 0.115. 



Experiment 


Central Value 


Error 


ALEPH jets 


0.125 


±0.005 


DELPHI jets 


0.113 


±0.007 


L3 jets 


0.125 


±0.009 


OPAL jets 


0.122 


±0.006 


OPALr 


0.123 


±0.007 


J/* 


0.108 


±0.005 


T 


0.109 


±0.005 


Deep Inelastic 


0.109 


±0.005 


Average 


0.116 


±0.005 



Table 1 - Experimental Values of a^{mz)- 



Input Values 


sin 2 9 w {m>z) 


m t 


(Xem(mz) 


m h 


li 


X 


-^-max 


0.2334 


95 GeV 


1/128.1 


ITeV 


ITeV 


193.2 


central sin 2 dw{ rn z)^t 


0.2327 


130 GeV 


1/128.1 


ITeV 


ITeV 


189.7 


central a em (mz) 


0.2334 


95 GeV 


1/127.9 


ITeV 


ITeV 


192.9 


-A- central 


0.2327 


130 GeV 


1/127.9 


ITeV 


ITeV 


189.4 


-^-min 


0.2320 


160 GeV 


1/127.7 


ITeV 


ITeV 


185.4 


X max : [i,m h < 500 GeV 


0.2334 


95 GeV 


1/128.1 


500 GeV 


500 GeV 


190.6 



Table 2 - Sensitivity of X to various inputs: 
/U, mh < 1 TeV unless otherwise stated. 
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